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We investigate models of (1 + d)-D Lorentzian semi-random lattices with one random 
(space-like) direction and d regular (time-like) ones. We prove a general inversion formula 
expressing the partition function of these models as the inverse of that of hard objects 
in d dimensions. This allows for an exact solution of a variety of new models including 
critical and multicritical generalized (1 + 1)-D Lorentzian surfaces, with fractal dimensions 
dp = fe+1, k = 1, 2, 3, as well as a new model of (l+2)-D critical tetrahedral complexes, 
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from this solution. We finally establish a general connection between (1 + d)-D Lorentzian 
lattices and directed-site lattice animals in (1 + d) dimensions. 
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1. Introduction 

The study of random lattices is an important subject with relations to many areas of 
physics such as quantum gravity or fluid membranes. One of the most efficient descriptions 
of such random lattices is through matrix models which generate discrete two-dimensional 
random surfaces in the form of tesselations made of tiles of arbitrary valences (for reviews 
see and references therein). The archetypical example is random triangulations pe- 
rused in the context of 2D quantum gravity. Analogous tesselations can be considered in 
higher dimension but unfortunately no such powerful tool as matrix models is available, 
making the subject quite difficult. Recently a new type of random lattices was introduced 
referred to as Lorentzian random lattices in which one particular (time) direction is 
regular while the other (space) ones are random. This allows to view a, say, (D + 1)- 
dimensional Lorentzian lattice as the time evolution of a D-dimensional usual random 
lattice. More generally, we may as well consider semi-random lattices, with a number, say, 
D of random directions and d of regular ones. We shall refer to these lattices as (D + d)-D 
Lorentzian lattices. These can be viewed as interpolations between regular lattices (D = 0) 
and fully random ones (d = 0). By a slight abuse of language, we shall refer to the d regular 
directions as "time directions" and to the D random ones as "space directions" , although 
our analysis is purely statistical (no dynamics). The cases considered so far correspond to 
d= 1 and D = 1 @[| or 2 [§. 

In the present work, we focus on the case of D = 1 and arbitrary d > 1, referred to 
as (1 + d)-D Lorentzian-type lattices in the following. For such models, we will derive a 
powerful inversion relation, expressing the partition function for such a (l+d)-D Lorentzian 
lattice as the inverse of the partition function of some hard object model on the time-like 
d-dimensional regular lattice. More precisely, this relation takes the general form 

ZiMU}) = wrm (L1) 

where Zi + d({ti}) is the partition function of our (1 + d)-D Lorentzian lattices with ac- 
tivities ti per tile of type i, and Z%({— ti}) denotes the partition function of hard objects 
on the time-like d-dimensional regular lattice with activities —t{ per hard object of corre- 
sponding type i. This relation ( |1 . 1| ) is very general, and holds for any dimension d and any 
fixed (i-dimensional (but not necessarily regular) time lattice, any type of tile i, and even 
with time-dependent activities £j. We may view the inversion relation ( |1 . 1| ) as some type 
of "boson-fermion" correspondence, relating the partition function of weakly interacting 
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bosons, namely the tiles of the lattice, to the inverse of that of locally interacting fermions, 
namely the hard objects with nearest neighbor exclusion. 

As a first application, we show how to relate the partition function of (1 + 1)-D 
Lorentzian triangulations to that of the hard dimer model on a line, which allows for a 
very simple solution of the problem in terms of a 2 x 2 transfer matrix, as opposed to the 
previous solutions relying on transfer matrices of infinite size fU Q || . This new approach 
allows to build and solve many more models of (1 + 1)-D Lorentzian surfaces, such as 
those made of larger (2(z + l)-gonal) tiles, in connection with hard multimers on a line. In 
particular, we are able to reach new multicritical points for these surfaces, displaying new 
large scale universal properties, and fractal dimension Dp = k + 1, k = 1, 2, 3, .... 

Going to higher dimension, we then introduce a model of (1 + 2)-D Lorentzian tetra- 
hedral complexes, i.e. semi-random lattices made of tetrahedra, and apply our inversion 
relation to express its partition function in terms of that of 2D hard hexagons solved by 
Baxter 0. As an outcome, we immediately obtain the large scale behavior of these new 
semi-random lattices. 

Analogous relations between (1 + d) -dimensional problems and <i-dimensional nearest- 
neighbor exclusion models have already been found in the context of directed-site lattice 
animal enumeration (DSAE) problems This suggests the existence of a connection 

between Lorentzian-type (1 + d)-D Lorentzian lattices and (1 + d) -dimensional DSAE. We 
will establish such a connection in which animals will appear as a particular subclass of 
Lorentzian lattices. Our inversion relation actually provides an alternative and more direct 
derivation of the equivalence between DSAE and nearest neighbor exclusion models. 

The paper is organized as follows. In Section 2, we discuss the case of Lorentzian 
(1 + 1)-D surfaces. We first derive the inversion formula by focusing on the simplest 
case of (1 + 1)-D Lorentzian triangulations (Section 2.1). This allows us to rederive very 
simply some of the known properties of these surfaces. In Section 2.2, we extend the 
inversion formula to include (1 + 1)-D Lorentzian surfaces made of time-like 2(i + l)-gons 
with activity £j i = 1, 2, 3..., now corresponding to hard (i + l)-mers on a line, and derive 
the corresponding thermodynamic partition function and loop-loop propagator. The first 
application concerns the case of surfaces made of only one type of such tiles (say 2(k + 1)- 
gons, with fixed k) and is discussed in Section 2.3. Next we show in Section 2.4 how to 
obtain multicritical models by fine-tuning the various activities ti. For these models, we 
compute the corresponding scaling exponents as well as universal scaling functions. Section 
3 is devoted to the study of (1 + 2)-D Lorentzian tetrahedral complexes. We first define 
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the model in Section 3.1 in terms of plaquettes living in tubes of hexagonal section. We 
then apply in Section 3.2 the inversion relation to obtain the critical behavior of our model 
in terms of that of hard hexagons at the Lee- Yang edge singularity point. In Section 4, 
we make the connection between our models and directed-site lattice animals, both in the 
(1 + 1)-D case (Section 4.1) and in the (1 + 2)-D one (Section 4.2). Finally, Section 5 is 
devoted to a discussion of the (infinite) transfer matrices for our models. We first derive 
in Section 5.1 the various (1 + 1)-D transfer matrices in terms of the corresponding finite 
ones for hard objects on a line. The most technical cases of this discussion are treated 
in Appendix A. We then use in Section 5.2 the equivalence to hard objects to construct 
more general parametric families of mutually commuting transfer matrices corresponding 
to integrable models containing our semi-random lattice models as particular points. We 
gather a few concluding remarks in Section 6. 

2. Critical and Multicritical Models of Lorentzian (1 + 1)-D Surfaces 

2.1. Inversion principle: (1 + 1)-D Lorentzian triangulations vs ID hard dinners 

In this section, we introduce the fundamental inversion formula relating the partition 
functions of Lorentzian-type semi-random lattices and hard objects in one less dimension. 
For simplicity, we specialize here to the simplest model of pure Lorentzian triangulations 
in (1 + 1)-D, which corresponds to hard dimers on a line. 




T 



space space 



Fig. 1: A typical (1 + 1)-D Lorentzian triangulation together with its dual, 
made of T time-slices. The triangulation is regular in the "time" direction 
and random in the "space" direction with an arbitrary succession of up and 
down triangles in each slice. Triangles of neighboring slices may be paired 
so as to form time-like lozenges, such as the shaded one in the figure. These 
elementary building blocks translate into vertical edges in the dual picture. 
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We start with the partition function Zt{£) of Lorentzian triangulations with T 
time slices and an activity t per pair of neighboring triangles sharing a space-like edge (see 
Fig.l). The corresponding "time-like" lozenges form the building blocks in the construction 
of the surfaces. In the dual picture, they translate into time-like (vertical) edges connecting 
successive space-like (horizontal) lines. 
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Fig. 2: A configuration of (1 + 1)-D Lorentzian triangulation in the dual 
picture, together with its left vertical projection. The latter is obtained by 
letting the foreground of the triangulation, made of the leftmost edges, slide 
horizontally all the way to the vertical line on the left. This projection clearly 
defines a hard dimer configuration on the vertical line. 

The idea behind the correspondence to hard dimers is to decompose the surface con- 
figurations according to their left vertical projection defined as follows. Let us allow for 
the vertical edges to slide along the horizontal direction without passing one-another both 
within the same time-slice and between two consecutive ones, thus preserving the relative 
positioning between them. We then single out those vertical edges which can be taken all 
the way to a vertical line on the left without moving the others. The resulting configura- 
tion of edges along this vertical line constitutes the left vertical projection of our surface. 
Clearly, it defines a hard dimer configuration of the vertical integral segment [0, T] C 2 
as depicted in Fig. 2. To get dimers we simply view each edge as linking its two endpoints. 
The hardness simply means that any integer point in [0, T] belongs to at most one dimer, 
translating into a mutual avoidance of dimers. This is due to the fact that on any two 
consecutive time-slices, no more than one edge can be projected. 
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We therefore write the partition function Zj-(t) as 



Z T {t) 



t 



\C\ y{C) 



(t) 



(2.1) 



hard dimer config. C 



where the sum extends over all hard dimer configurations C on the integer segment [0, T] 

(C) 

(including the empty one), and where Zj, (t) is the restricted partition function involving 
those configurations having projection C, from which we have factored out the weight 
of the projected part, \C\ denoting the number of dimers in C. More generally we have 
the relations 



T 

D 



(a) 



(b) 



(c) 



Fig. 3: By completing an arbitrary Lorentzian triangulation with a hard 
dimer configuration D (a), we build a larger triangulation (b) whose projec- 
tion C contains D (c). With this procedure (with fixed D) we build exactly 
once all Lorentzian triangulations whose projection C contains D. 



t \ D \z T (t)= tlCl 4 C) (t) (2-2) 

CDD 

valid for any hard dimer configuration D (eqn. (|2~T|) corresponding to D = 0). This ex- 
presses the fact that by completing any configuration of Lorentzian triangulation with T 
time slices (in the dual picture) by a given left column of vertical edges (corresponding to 
a hard dimer configuration D), one builds exactly once each configuration having a pro- 
jection containing D, i.e. having D as a sub-configuration (see Fig. 3). This latter relation 
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is easily inverted using the celebrated Mobius inversion! formula, leading to 

4 C) (t) = (-l) ]Dl - lC] t^Z T (t) (2.3) 

DDC 

Noting that Zx(t) factors out of the sum on the rhs, we finally get 

DDC 

Picking C = 0, we arrive at our fundamental inversion relation 



1 

Zff{-t) 



where 

z^ d ( z ) = J2 zlDl ( 2 - 6 ) 

hard dimer config. D 

denotes the standard partition function for hard dimers with fugacity z per dimer. As we 
already mentioned in the introduction, this relation is a generalization of the boson-fermion 
correspondence relating for instance the partition function 1/(1 — t) of free bosons with 
fugacity t per particle on a point to the inverse of 1 + z = 1 — t of that of a fermion with 
fugacity z = — t on a point. 

The formula (|2.4j) also implies, upon substituting (|2.6[) that 



where 



Z[O d T]\c( 

Z^ d {-t) 



ft,\cW=EW |DHC| (2-8) 

DDC 

is the hard dimer partition function on the segment [0, T] minus the occupied points of C. 

The above construction is very general. In particular, our inversion formula ( |2.5| ) 
holds also if we attach an extra fugacity w s per edge inside the time-slice s resulting in a 
total weight z s = —tw s per dimer in the slice s in the hard dimer language, allowing in 
particular to recover the loop-loop propagator (partition function with a fixed number of 



3 Note that the Mobius inversion involves usually the opposite order relation C but it works 
similarly for the order relation D used here. 
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triangles at times and T) , and correlation functions for the numbers of triangles in given 
slices. It will also be extended in Section 2.2 below to the case of surfaces made of larger 
building blocks corresponding to hard multimers on a line. Finally in Section 3 we will 
extend it to higher-dimensional semi-random lattices. 

With this construction we have reduced the (1 + 1)-D Lorentzian gravity-type problems 
to that, much simpler, of hard objects on a line. 

To complete this section, let us re-derive the partition function of pure Lorentzian 
triangulations Zx{t) from the hard dimer equivalence. The partition function Z!^ d (z) is 
easily computed by use of a 2 x 2 transfer matrix T between successive segments in either 
empty or occupied states: 

Zt*{z)=v*f r *v = {l 0)Q ^^(J) (2-9) 

We may read directly from this the asymptotic large T behavior of Zx{t). Indeed, ( |2.9j ) 
is dominated for large T by the largest eigenvalue A_|_ of T, satisfying of A 2 — A — z = 
(with solutions X±(z) = (1 ± y/1 + Az)/2). Therefore we find that for large T 

Zr{t) = z J { _ f) ~ v T (t) where fj,{t) = (2.10) 
is the smallest solution of 

H = l + tn 2 (2.11) 

namely the celebrated generating function of the Catalan numbers fi(t) = (1 — y/l — At) / (2t). 
More precisely, we have 

Z^ d (z) -- 

leading to 

z T (t) 

where we have set q = \ftji{t) (or equivalently 1/y/i = q + 1/q). Expanding the de- 
nominator in ( |2.13| ) we immediately read the eigenvalues of the infinite transfer matrix 
for Lorentzian triangulations (q + ^)q 2n+1 , n = 0,1,2,..., thus recovering the harmonic 
oscillator eigenvalues 0] . 

As a final remark, we may interpret fj,(t) as the partition function of a sub-class of 
triangulations on the semi-infinite time interval [0, oo). To get a finite partition function it 
is sufficient to demand that the projection of the configurations is either the vacuum or the 



A + (z) T+2 -A_(z) r + 2 
X+(z)-X-(z) 



(2.12) 



2\T+2 



1-g 2 + g 
1 + q 2 1 — g 2(T+2) 



(2.13) 
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Fig. 4: A typical Lorentzian triangulation with projection Co, made of a 
single dimer in the lowest position. This induces a staircase-type boundary 
condition on the left side of the triangulation, namely that any edge in a given 
time-slice must have at least one edge on its left in the slice just below. 

hard dimer configuration Co made of a single dimer in the first slice (at time 0, see Fig. 4). 
This implies a so-called "staircase" boundary condition on the left of the configurations, 
namely that an edge occurs in the slice s only if an edge already occured on its left in the 
slice s — 1 below. Indeed, we write this partition function as 



As explained in Ref. 0] and depicted in Fig. 5, the configurations above are in one-to- 
one correspondence with discrete Random Walks with steps of ±1 on the integral half-line 
[0, oo), starting and ending at the origin, and with a weight t per ascending step. 

As we will see in the following sections, all the above equivalences will nicely generalize 
to more involved cases. 

2.2. Generalized (1 + 1)-D Lorentzian surfaces vs ID hard multimers 

In this section, we introduce generalized discrete Lorentzian surfaces made out of 
various tiles, including lozenges, hexagons, octagons, etc... To be more precise, we wish 
to compute the partition function ZT{{ti}) of surfaces which in the dual picture look like 
Fig. 6, with vertical edges of arbitrary length with a fugacity tj per edge of total length i, 




(2.14) 



= l + t lim - 



T 
,T-2 



l+tfl 2 = n 



Fig. 5: The one-to-one correspondence between "left staircase boundary" 
triangulations and random walks from the origin to itself on the integer half- 
line. We first rewrite the configuration of edges forming the triangulation as a 
tree, by connecting each edge to that sitting just below on its left (connections 
are represented by thick horizontal lines here). By following the contour of 
the tree from the lower left branch to the lower right, the sequence of ascents 
and descents gives rise to the directed walk represented in thin solid line (to 
get a nice directed walk, we must first place each vertical edge at a horizontal 
position equal to total number of ascents and descents along the tree preced- 
ing it). Conversely, each such walk gives rise to a unique tree, therefore to a 
unique triangulation with the staircase boundary condition. This correspon- 
dence includes the empty triangulation with projection 0, corresponding to 
the walk of step. 

i = 1,2, 3... and taken on a time segment [0, T]. These edges are dual to particular time- 
like 2(i + l)-gons (see Fig. 6 (a)). The inversion formula (|2.5|) generalizes straightforwardly 
to this case. Firstly we now have the generalization of eqn. ( |2.2| ) 

\{t? {D) z T {{ t% }) = y: nc (c) zFau}) (2.i5) 

i CDD i 

where D (resp. C) denote hard multimer configurations on the integral segment [0, T] 
with rii(D) (resp. n^(C)) (z + l)-mers. Again in such a configuration hard multimers avoid 
one-another in that a given point may belong to at most one multimer. By inclusion of 
configurations C D D we mean that C contains all the multimers of D plus possibly others. 
Eqns. (|2.15| ) are easily inverted by the Mobius inversion formula to finally yield: 

1 



Zt{{U}) = 7TT (2-16) 



where Z^ m ({zi}) denotes the partition function of hard multimers on the integral segment 
[0, T] with a fugacity z% per {% + l)-mer. 
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(a) 



(b) 



Fig. 6: A typical generalized (1 + 1)-D discrete Lorentzian surface (a) made 
of time-like lozenges, hexagons, octagons, etc... respectively weighted by 
tii *2, £3, ••• Its dual (b) is made of vertical edges of length 1, 2, 3, ... extending 
over several time-slices. The corresponding left vertical projection is nothing 
but a configuration of the hard multimer model on the integer segment [0, T]. 

Again, in the hard multimer language we simply have to diagonalize the corresponding 
transfer matrix. Let us fix for definiteness a maximal length k for the edges, corresponding 
to a model of di-, tri-, (k + l)-mers. This truncates the transfer matrix to a size 
(k + 1) x (k + 1), where the (k+ 1) possible states correspond to the empty state, the state 
occupied by the lowest monomer of a (i > l)-mer, the state occupied by the second lowest 
monomer of a (i > 2)-mer, ... The transfer matrix then reads 



/I 


1 


1 



















£3 
Z2 



\ 

and the partition function is expressed as 



1 






Zk- 



1\ 







(2.17) 



Z^({z l })=v t (T k ) T+1 v 
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(2.18) 



where v* = (1, 0, 0, 0). Again the large T asymptotics of ( 2.16|) are governed by the 
largest eigenvalue \({zi}) of T%, leading to 



Zt({U}) ~ V T ({U}) (2.19) 

where n({ti}) = 1/X({zi = —W}) is the smallest (in module) solution of the characteristic 
equation 

k 

n = i+j2 t ^ i+1 ( 2 - 20 ) 

Here again, \x may be interpreted as the partition function for the particular Lorentzian 
surfaces with semi-infinite time interval [0, oo), having a projection either empty (contri- 
bution of 1 in the rhs of fl2.20p ), or equal to a single multimer of either length i, extending 
over the i first time slices, i = 1, 2, k (contributions of tifi t+1 in the rhs of (|2.20|) ). 




Fig. 7: The one-to-one correspondence between "left staircase boundary" 
generalized discrete Lorentzian surfaces and random walks from the origin to 
itself on the integer half-line, with ascending steps of arbitrary length and 
descending steps of —1. As in the triangulation case, we associate to each 
"left staircase boundary" surface configuration a tree connecting each edges 
to that sitting just below it on its left. The walk is then defined as a walk 
along the tree, with the convention that it makes an ascending step of +i 
when going up along an {i + l)-mer, whereas it only makes descending steps 
of —1 (again to best represent it as a directed walk, we must place the edges 
at horizontal positions equal to the total number of ascents and descents 
preceding them). Conversely any such walk gives rise to a unique tree with 
branches of arbitrary integer lengths and therefore a unique "left staircase 
boundary" surface. 
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As before, these configurations are in one-to-one correspondence with discrete Random 
Walks from to on the integral half-line [0, oo), with possible ascending steps of +i, 
i = 1, 2, /c, weighted by ti, and descending steps —1 (see Fig. 7). The partition function 
[i of these walks clearly satisfies the self-consistent equation ( |2.20|) corresponding to a 
decomposition of the walks according to the length i of their first ascending step, with 

being the partition function of the rest of the walk going from position i back to the 
origin on the half-line. 

Another quantity of interest is the so-called loop-loop propagator Zx(i, j\{t m }) de- 
fined as the partition function of surfaces with T time slices and with fixed numbers 
i (resp. j) of tiles originating (resp. terminating) in the first (resp. last) slice. The 
inversion formula (|2.16| ) extends straightforwardly to include an extra weight w s at- 



tached to each portion of tile visiting the slice s = 0,1, ...,T. Choosing wo = x, 
wt = y and w 8 = 1 for s = 1,2,...,T — 1, we immediately get the generating func- 
tion Gt(x, y\{t m }) = J2i j>o x% V^T{i) j\{tm}) as the inverse of the partition function 
Gy m (x, y\{z m }) of hard multimers with weights z m = —t m for (m + l)-mers and extra 
weights x (resp. y) for the multimers originating (resp. terminating) at point (resp. T). 
To compute G^ m (x, y\{z m }) we make the simple observation that it is at most linear in 
the variables x and y as at most one multimer can touch the point (resp. T), and that 
it is symmetric in x and y. Writing 

G^ m (x, y\{z m }) = a-b(x + y) + cxy (2.21) 

we easily compute the coefficients a, 6, c by taking particular values of x and y. For x = 
y = 0, we forbid multimers to touch the extremities of the segment [0, T] thus effectively 
reducing it to [1,T — 1] and resulting in a = ZjP^ 2 ({ z m}) ■ For x = 1 and y = we only 
forbid the T point, with the result a — b = Zl^ 1 {{z rn \) . Finally, for x = y = 1 we simply 
have a — 2b + c = Z^ m ({z m }). This results in the following expression for the loop-loop 
propagator 

G T (x,y\{tm}) = — t ; — -r j—k 7—l u E—r ( 2 -22) 

i\ ,y\\ mS! z^ 2 -(x + y)(Z^ 2 -Z^ 1 )+xy(Z^ 2 -2Z^ 1 + Z^ m ) v ; 

where the Z^ 1 are all taken at the values z m = —t m . Explicitly expanding this as a series 
in x and y, we obtain 



i / zh™ \ x ^ fi\ fj\ ( (z^) 2 - z£ m z£™ 2 

z T ( i ,Mt m} )= w - i-^i Egg ( z- i (2 - 23) 



T-2 
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2.3. First application: hard (k + l)-mers and Fuss-Catalan numbers 

As a first simple application of our construction, let us consider surfaces made only 
of time-like 2{k + l)-gons, i.e. with duals made only of vertical edges of fixed length k. 
This corresponds to specializing the above to tj = S^t, where t is the weight per tile. 
In the limit of large T, this yields the thermodynamic partition function Z^\t) ~ fx 7 \t) 
where n(t) is the smallest solution (in module) of tfi k+1 = \i — 1, known as the generating 
function of the Fuss-Catalan numbers cl [IT 



Mt) = E c ? >r where c« = {|±iM; (,24, 

n>0 V ' 

( k) 

From the use of Stirling's formula, we get the large n behavior of these numbers: c n ~ 
[ik + l) fc+1 //c fc ) n /n 3 / 2 , which allows to show that for all k > 1, the function fi(t) displays 
a square root singularity /i| s i ng ~ \Jt c — t when t approaches the critical value 

*° = (^TT (2-25) 

Therefore the scaling limit of these models (for all k > 1) lies in the same universality class 
as that of pure Lorentzian triangulations, corresponding to k = 1. 

Note finally that according to the above equivalence with Random Walks, the function 
fj,(t) also generates the numbers of directed Random Walks on the integer half-line [0, oo) 
starting and ending at and with ascending steps of +k only and descending steps of —1, 
and with a weight t per ascending step. 

2-4- Second application: multicritical models of (1 + 1)-D Lorentzian surfaces and Patalan 
numbers 

In this section we show how to go beyond the generic square root singularity of pure 
Lorentzian surfaces and get more interesting critical behaviors. As usual, this can be done 
by fine-tuning the weights U in order to reach multicritical points. Indeed, we may reach 
a multicritical point of order (k + 1) by retaining di-, tri-, ... , (k + l)-mers and fine-tuning 
the activities ti, i = 1, 2, k in order for eqn. ( |2.20| ) to take the form 

2 )W .. = 1 _2(* + l) t 
v k ' k 
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namely by picking 



k + l\ ( It 



1, 2, ...k (2.27) 




all expressed in terms of the activity t = t± per lozenge. The values of the coefficients in 
( p.26|) are entirely fixed by the constant, linear and quadratic terms in ( |2.20| ), with relative 
values 1, —1 and t. This yields the thermodynamic partition function per time slice 

1-2^)*) (2.28) 

For instance, for k = 2, we reach a tricritical point by taking activities t per lozenge 
and — 1 2 /3 per hexagon, leading to fx = (1 — (1 — 3t)^)/t. The need for both positive 
and negative activities to reach a multicritical point parallels the case of ordinary random 
surfaces as solved by means of the one-matrix models, whose potentials display the same 
pattern of alternating signs for the activities per tiles. Note also that similarly to that case, 
we need to consider at least {k + l)-mers to reach a multicritical point of order (k + 1). In 
that respect, eqn.( |2.26| ) is the minimal realization of the (k + l)-critical point in which the 
requirement of multicriticality fixes all the t^s in terms of t = t±. The same point could 
be attained from any model involving also larger multimers leading to a characteristic 
equation of degree larger than k + 1. The corresponding (k + 1) -critical point involves then 
k — 1 relations between ti, tz, .... 

Up to the change of variables t = k(k + l)x/2, the function fi = ^2 n>0 Pn x n ( |2.28| ) 
is known as the generating function of so-called Patalan numbers [|TIJ 

which are all positive integers. The integrality of these numbers is clear from eqn. ( |2.27|) , 
since all the ti are integer multiples of x l . On the other hand, the positivity appears 
as non-trivial, as the ti have alternating signs. This strongly suggests a possible purely 
combinatorial reinterpretation of our multicritical partition functions. 

The singularity of ( |2.28 ) leads to a new scaling behavior for the corresponding dec- 



orated surfaces, when t approaches the multicritical point t c = k/(2(k + 1)). Indeed, to 
get a proper scaling limit of the multicritical partition function Z^\t) — Zx({ti}) ~ /U T , 
corresponding to the fine-tuning ( |2.27| ) and with \i as in ( |2.28| ), we must set 

T=- t = t c (l-a k+1 A) (2.30) 
a 
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where r is the renormalized time lapse and where A is the renormalized activity per unit 
area of tile ( "cosmological constant"), with a — > in the scaling regime. In this regime, 
(/i(t)//i(t c )) T — ► e ~ TAk+1 and more generally Z^\t)/Z^\t c ) becomes a universal scaling 
function of rA 7 =+ T which we will determine below!. Introducing the area A = Yl^i^i 
where iVj is the total number of 2{i + l)-gonal tiles forming the surface, we see that the 
total Boltzmann weight is proportional to t A . Writing t A ~ t A e~ A r ° when t — > t c , and 
substituting Q2.3Q ), we see that the parameter A is conjugate to the renormalized area of 



the surface, defined as A = a h+1 A. We immediately deduce that A ~ T h+1 and get the 
scaling behavior of the area of the surface 

A~T k+1 (2.31) 

in terms of the time lapse T when T is large. This determines the "fractal dimension of 
space-time" to be cLf = k + 1. In the directed Random Walk picture of previous section, 
this allows to obtain the scaling behavior of the vertical extent or gyration radius R = T 
of the walks on [0, oo), in terms of their total length L = 2A, namely R ~ L v with the 
exponent 

z,= -L = — L_ (2.32) 
d F k + 1 K 1 

This is also the exponent of the correlation length £ ~ (t c — t)~ u appearing for instance 
in the correlation of the numbers of tiles A/"(Ti), MiT^) in two given time slices at times 
Ti,T2, namely (Af(Ti)Af(T2)) ~ e « . The partition function of the walks of fixed 
length L is identified with the coefficient of in the expansion of n(t) ( |2.28| ), and behaves 
as (t c )~^L a_3 with the susceptibility exponent 

2A; + 1 . 
a = (2.33) 

Equivalently, in the language of random surfaces, the coefficient of t A in the expansion of 
li(t) ( [2.28 ) represents the thermodynamic partition function Za of fine-tuned semi-random 



surfaces of fixed area A and with fixed projection either empty or equal to Co, namely 
reduced to either zero or a single dimer. Therefore the exponent a ( |2.33|) is interpreted as 
a configuration exponent for these objects, namely Za ~ (t c )~ A A a ~ 3 for large A. Note 



4 Note that this scaling function involves /x and also the other roots of the polynomial equation 
( 2.26| ) as they all merge at the multicritical point. 
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finally that, as expected, the exponents a and v above obey the hyperscaling relation 
2 — a = dv in d = 1 dimension. 

To conclude this section, let us derive the explicit form of the finite-time partition 
function zip (t) for the fine-tuned semi-random surfaces and its scaling limit, together 
with that of the corresponding loop-loop propagator. Starting from the inversion formula 
Q2.16|) , we are left with the calculation of the partition function of fine-tuned hard multimers 

ryhmlk] -i 

Z T expressed as 

Z hm[k] =v t^T+l v (234) 

where the (k + 1) x (k + 1) transfer matrix reads as in fl2.17f) with z% = —t\ given by ( |2.27|) , 
and where v l = (1, 0, 0, 0). The eigenvalues Xj, j = 0, 1, k of this transfer matrix are 
simply the inverses of the solutions /ij to the multicritical characteristic equation ( |2.26|) , 
reading 

N = -(1-^A) with^ = e 2i ^r, A= (1- 2(/C + 1) t)^r (2.35) 
2t k 

The partition function ( |2.34j ) is expressed as a linear combination Z^ m ^ = J2o<j<k a jl (f x j) T+1 i 
where the coefficients aj may be obtained by explicit diagonalization of 7^. After some 
algebra we find the rather simple expression 



which reduces to ( |2.12| ) for k = 1. Note that it satisfies the recursion relation Z^j fc ' = 
2?hrn[k] _ j^^ =1 tiZ T ™\ k ^ with the fine-tuned values of ti ( |2.27| ). The solution ( |2.36j ) is the 
unique function obeying this recursion relation with the initial conditions: = 0, 

i = 2, 3, k + 1 and Z^ _ ^ p or jn us tration, let us display the first few values of 

Zhmlk] r i • rr-i 

T 1 ' tor non-negative 1 : 

ryhm\k\ -. 

z =1 

^hm[k] ^ 

z hm[k] = 1 _2t+ 2 ( k ~ 1 ) f 2 ( 2 - 37 ) 

2 3k 



1 (It 



T+2 k 



7 hm 



Z'z m = 1 - 3t + - t z - 



7k -4 2 (fc-l)(fc-2). 



3k 3k 2 
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By use of the inversion formula ( 2.16 ), the expression (|2.36| ) finally leads to the partition 



function of multicritical semi-random surfaces 



4%)=(^) ^ k ~ uj (2-38) 

'j=0 (l-^A) T + 2 



k 



T+2 



2t J V fe ^~ 

2^j=o (i- 



Note that as t — > t c = k/ (2(k + 1)), A — > 0, the partition function function ( |2.38j ) tends to 
a finite limit, as the denominator is of order A fc . This yields a finite ratio 

Z T ] ( t )_ n i\ ft c \ T+2 fT + l + k 



(*+!)( t) ( V ) ^. (2-39) 

(1-^'A) T + 2 



We may now perform the scaling limit ( |2.30| ) of this ratio, that tends to a universal scaling 
function of the variable x = rA^+i, namely 

lim £ W = „ , . , , x = tK—^ 2.40 

where we have used the generalized hyperbolic function h\ of order k + 1, member of the 
family h p (x), p = 0, 1, 2..., fc, defined by 

fc 

hp(a;) = ^u; p V ,iie (2.41) 



with cj as in ( 2.35 ). The functions h p are cousins of the Mittag-Leffler function 



generalizing the hyperbolic sine and cosine (case k = 1, p = 1, 2 respectively). For k = 1 
indeed, with = 2 sinh(x), ( |2.4U| ) allows to recover the usual scaling function xj sinh(x) 

i- 

We can fix the normalization of the partition function in order for its scaling limit to 

read 

Zp (A) = lim - T \V = V ; t 2.42 

1 } (k + 1)! a^o (k + 1) T klhiirA^) 

and accordingly the critical value z\ (0) = l/r k . In particular, this allows to recover for 
k = 1 the partition function z| 1] (A) = \/A/ sinh(rv / A) of g]. 

Let us now turn to the loop-loop propagator (|2.23|) . To define a sensible scaling limit 
of this quantity we need to take 

i = ^ 3 = ^ (2-43) 
a a 
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and to substitute in (|2.23| ) the value of the hard multimer partition function 

Z$™ a cx —±— hl (TA&) (2.44) 

with hi as in ( |2.41| ). Indeed, except for the prefactor 1/Zj,™ 2 which we normalize as in 
( |2.42|) , only ratios of Zy m, s enter the expression ( |2.23| ) so we may forget about all the 
prefactors independent of T as indicated in ( 2.44j) by the proportionality symbol. Eqn. 
(|2.23|) finally takes the scaling form 



G^UM) = 4"(A) /„ ( 2^L lia ^lM^l*i| (2 . 46 ) 

with Z^\A) as in ( |2.42| ), with Io(2x) = J2 r >o x 2r /(rl) 2 the modified Bessel function, and 
where the hyperbolic functions h p = ^(rA^+i) are defined in ( p.41|) . For completeness, 
let us display the hyperbolic functions entering the loop-loop propagator for the first few 
values of k = 1, 2, 3. 

k = 1 : hi(x) = hs(x) = 2 sinh(:c) h<z{x) = 2cosh(x) 

k — 2 : h m (x) = e x + 2e~ ? cos( — x + 2m— ) m = 1, 2, 3 
fc = 3 : /ii(x) = 2(sinh(x) — sin(x)) h,2(x) = 2(cosh(x) — cos(x)) 
h%(x) = 2(sinh(x) +sin(x)) 

With these, we recover in the k = 1 case the results of H. 

Let us finally comment on the form of the rescaled partition function z\ k \A) ( |2.42|) . 

In terms of the rescaled variable x = tA 7 h 7T , it has poles situated at the complex non- 

i 

vanishing zeros of h±(x), all of the form Xj >m = —u^am 1 , J ' = 0, 1, 2..., k and m = 1, 2, 3, ... 
with a m real positive. This allows to write the partition function as 

Z[ k] (A) = ri ^ r - (2.47) 

r fe+1 n~=i(l + ^) 

When k = 1 we simply have a m = 7r 2 m 2 , m = 1, 2, 3.... For arbitrary k there is no such 
simple form, but asymptotically one can show that 

a m = — — z [m+— - + 0(e" em ) (2.48) 
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3. (1 + 2)-D Lorentzian tetrahedral complexes and the hard hexagon model 

In this section, we consider three-dimensional random objects made of tetrahedra, 
corresponding to (1 + 2)-D Lorentzian tetrahedral complexes, with two (time-like) regu- 
lar directions and one (space-like) random one. This is to be contrasted with the more 
complicated situation of discrete (2 + 1)-D Lorentzian manifolds used in the context of 
3-D gravity |J, which have one (time-like) regular direction and two (space-like) random 
ones. Still, by slightly breaking the symmetry in the two time directions, we shall be able 
to view the configurations of our model as the time-evolution of some particular random 
triangulations, namely the Lorentzian ones. 

3.1. The plaquette model for Lorentzian tetrahedral complexes 



Fig. 8: A sample hexagonal plaquette configuration (a). There are arbi- 
trarily many plaquettes in each tube. Neighboring plaquettes cannot cross 
each other, neither within the same tube, nor between adjacent tubes. Such 
configurations are dual to tetrahedral complexes, made of tetrahedra filling 
triangular tubes dual to the hexagonal ones. We have represented an elemen- 
tary diamond-shaped building block (b) dual to a plaquette of (a), and made 
of six tetrahedra sharing an edge. 

If we view the pure (1 + 1)-D Lorentzian triangulations in the dual picture as a 
regular array of slices in which edges may freely slide horizontally provided they do not 
cross one-another, a natural (1 + 2)-D generalization consists in having a regular 2D array 
of horizontal tubes in which plaquettes may freely slide horizontally, provided they do 
not cross one-another. More precisely, we consider here the situation depicted in Fig. 8 




(a) 
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(a), where the tubes have a regular hexagonal section, and their array forms the two- 
dimensional hexagonal lattice. The plaquettes are filled hexagons orthogonal to the tubes 
and with the same section. A configuration is characterized by an arbitrary arrangement 
of such plaquettes in the tubes, where we take into account the relative ordering of the 
plaquettes with respect to one-another. In practice, a given plaquette only sees its nearest 
neighbors within the same tube and those in the six adjacent tubes, through their edges. 
For definiteness we attach a weight t per plaquette and consider only a finite array of tubes 
with transverse size T\ XT2 and with a bulky shape. With this definition, we generate three- 
dimensional objects with two (time-like) regular directions and one (space-like) random 
one. 

Like in the case of triangulations where the edges are dual to pairs of triangles form- 
ing lozenges, the above hexagonal plaquettes are dual to a diamond-shaped dodecahedral 
volume made of 6 tetrahedra glued around a common edge, as shown in Fig. 8 (b). The 
role formerly played by triangles is now played by tetrahedra: these live inside tubes with 
triangular section (dual to the original hexagonal tubes). Each such triangular tube is tes- 
sellated by tetrahedra of three possible kinds: each tetrahedron has 1 vertex on each edge 
of the tube, and the three kinds correspond to the three possible positions of the fourth 
vertex on one of the three edges of the tube§. Therefore exactly one edge of the tube 
contains one (space-like) edge of the tetrahedron. This edge is also common to the five 
other tetrahedra which together with the above form a diamond-shaped building block. 
All tetrahedra in a tube are glued along their triangular faces to two others in the same 
tube and two others in the two adjacent tubes sharing the abovementioned space-like edge. 

It is interesting to view our tetrahedral complexes as a random version of some regular 
lattice, in the same way as the Lorentzian triangulations are deformations of the triangular 
lattice. To construct a regular lattice, we simply need to fill the triangular tubes with 
tetrahedra in a regular way, say by taking a succession of tetrahedra of the first, second, 
third, first, etc... kinds along each tube. The other choice of chirality (first, third, second, 
first, etc..) is completely equivalent. The resulting lattice is nothing but a FCC lattice (see 
Fig. 9). Indeed, the latter can be viewed as a regular arrangement of octahedra completed 
by tetrahedra. Each octahedron can be decomposed into four tetrahedra by adding a 
diagonal edge (see Fig. 9). One possible elementary cell is obtained by first considering 

5 Note that this generalizes straightforwardly the (1 + 1)-D situation, where the triangles in 
each slice are of two kinds: those pointing up and those pointing down. 
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Fig. 9: The regular tetrahedral complex is nothing but the FCC lattice, made 
of octahedra supplemented by tetrahedra as shown. To get only tetrahedra, 
we simply add a diagonal edge inside all octahedra as indicated. This allows 
to decompose each octahedron into four tetrahedra. The building block of our 
model is composed of 6 tetrahedra, 2 from the original FCC lattice, sharing 
one edge, and 2 from each decomposed adjacent octahedron. 

two original tetrahedra of the FCC lattice sharing one edge, and then by completing 
them by the four adjacent tetrahedra (sharing the same edge), two from each neighboring 
octahedron. These 6 tetrahedra form the diamond-shaped building block of our model. 
The two choices of chirality correspond to the two possible choices of diagonal edge added 
inside the octahedra. In the dual plaquette language, the FCC lattice corresponds to 
the following regular arrangement of plaquettes: we decompose the (triangular) lattice 
of hexagonal tubes into three sublattices denoted, say 1, 2, 3, such that no two adjacent 
tubes belong to the same sublattice. Plaquettes are then arranged in equidistant successive 
parallel planes in such a way that plaquettes in successive planes cover the three sublattices 
1, 2, 3, 1, ... alternatively. 

The restriction of any given tetrahedral complex to a plane generated by the lon- 
gitudinal direction of the tubes and one of the three directions of the triangular lattice 
section of the tubes is nothing but a particular (1 + 1)-D Lorentzian triangulation (see 
Fig. 10). We may therefore view the configurations of our model as the time evolution 
of such triangulations in successive parallel planes (with a total of T 2 plane slices), each 
of which extends over a fixed time lapse T\ as shown in Fig. 10. The constant time-lines 
for the Lorentzian triangulations in two successive planes are shifted by half a period of 
the triangular lattice of tubes: indeed, these are nothing but edges of a same triangular 
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Fig. 10: The horizontal section of a sample (1 + 2)-D tetrahedral com- 
plex with vertical size T2 and horizontal size T\ is a (1 + 1)-D Lorentzian 
triangulation over a time lapse T\. 

tube. Each triangle in a given plane is linked to a vertex of the triangulation of the next 
plane, belonging to the same tube. These "face-to-vertex" tetrahedra clearly coexist with 
"vertex-to-face" ones, but also with "edge-to-edge" ones, namely with exactly one edge in 
each plane (one space-like and one time-like), and only time-like triangular faces. 
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Fig. 11: The dual configurations of two consecutive plane sections of a 
sample (1 + 2)-D Lorentzian tetrahedral complex. The lower (resp. upper) 
configuration is represented by dashed (resp. solid) lines. The constant time 
lines of the two configurations are shifted by half a period. We have encircled 
the three types of vertices of the resulting mixed graph, in one-to-one corre- 
spondence with the three types of tetrahedra filling the slice between the two 
planes (see text). 

The dual picture allows for a better understanding of the tetrahedral complex between 
two consecutive planes ||. Indeed, the dual configurations of the triangulations in the lower 
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and upper planes may be represented as in Fig. 11 in dashed and solid lines on the same 
picture, with their time-lines shifted by a half-period. The three types of vertices (solid 
trivalent, dashed trivalent and mixed tetravalent) of the resulting mixed graph are in one- 
to-one correspondence with the three types of tetrahedra (face-to-vertex, vertex-to-face 
and edge-to-edge) in the slice between the two planes. 

3.2. Equivalence with hard hexagons and critical behavior 

We now wish to evaluate the partition function Zx>{t) for the Lorentzian tetrahedral 
complexes in the plaquette formulation inside hexagonal tubes with a total section made of 
a compact connected domain T> of the hexagonal lattice. For each plaquette configuration, 
we define a left projection obtained by letting the leftmost plaquettes slide along the tubes 
all the way to their left end. The projections thus obtained are in one-to-one correspon- 
dence with configurations of hard hexagons on the two-dimensional domain T>. By hard 
hexagons, we mean that if an hexagon of T> is occupied, then its immediate neighbors must 
be emptyi. Repeating the arguments of Section 2.1, we arrive at the following inversion 
formula: 

where we have denoted by Z^ h (z) the partition function of the hard hexagon model on 
the domain T> with an activity z per occupied hexagon. We also have the more general 
formulas generalizing ( |2 . 7| ) : 



^'(t) = ^4 (3.2) 



where Z-p (t) is the partition function for the configurations with left projection C C T> 
and with the weight of the projection removed, while Z^ c (z) is the hard hexagon partition 
function on the domain T>\C, i.e. the set of hexagonal tubes in T> with no common edge 
with hexagons of C. The latter quantities ( |3.2| ) have a well-defined thermodynamic limit, 
obtained by keeping C fixed, and letting the domain T> become infinite. In particular, if 
C is reduced to a single plaquette Co, the quantity z!p°\t) tends to the thermodynamic 



6 In the standard picture, the hardness constraint is implemented by a weaker no-overlap con- 
straint for slightly larger hexagons, obtained from our hexagonal plaquettes by rotating them by 
7r/6 and dilating them by v3- These larger hexagons are also the projections of our diamond- 
shaped building blocks. The two ways of implementing the hardness constraint are clearly 
equivalent. 
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density of occupied hexagons Phh(z = —t) of the hard hexagon model. Note that this 
partition function Zj^ (£) is the generating function of the number Ny of tetrahedral 
complexes with V plaquettes (i.e. tetrahedral complexes with volume V if we assign a unit 
volume to the diamond-shaped building block) and projection a single one Cq. Note that 
the constraint that the projection be reduced to a single plaquette amounts to generalized 
"staircase" boundary conditions. 

We may now use Baxter's exact solution of the hard hexagon model [[7] to obtain the 
critical behavior of Zj^ (t). The hard hexagon is known to have two critical singularities 
at a positive and a negative values of z. The singularity of our problem occurs at the 
negative value z c = ( 1 ~ 2 V ^ ) 5 = 11 ~ 2 5v ^' of z (corresponding to a positive t c = —z c ), where 
the model is known to belong to the universality class of the Lee- Yang edge singularity 
|T3|1 [|14]], itself described by a conformal field theory of central charge c = —22/5 [[IS"] . The 



thermodynamic partition function per hexagon Khh(z) behaves 

K hh (z) = k + m(z - z c )i + 0(z - z c ) (3.3) 



corresponding to the critical exponent a = 2 — 5/6 = 7/6, and we also have 

Phh(z) = ■^-LogKhhi?) ~ \ — (z ~ z c )~* (3.4) 

dZ D Kq 

From this result, we immediately deduce that Zrp (t) ~ (t c — t) -1 / 6 , therefore the number 
Ny of tetrahedral complexes with V plaquettes and projection a single one behaves for 
large V as 

N v ~ const. % ' (3.5) 

V 'e 

Returning to the case of a finite but large domain D with area T\T%, we may write the 
partition function Zj)(t) at leading order in T1T2 as 



ZAi) - ^fw^ ~ KflT ' x (1 + > - i)lrT " (3 ' 6) 



Letting the time lapses Tj scale like Tj = — , a a small parameter, we see that we must 
approach the critical point as 

t = t c (l-a^A) (3.7) 

where A is the renormalized fugacity per tetrahedron, conjugated to the renormalized 
volume V = Va 12 ^ 5 . We deduce that for large Tj 

V^(TxT 2 )l (3.8) 
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hence the fractal dimension of the (1 + 2)-D Lorentzian tetrahedral complexes is dp = 
12/5. In other words, the transverse area T1T2 explored by the tetrahedral complexes 
with V plaquettes scales like T1T2 ~ V®, to be compared with the vertical length T of 
triangulations with A lozenges scaling like T ~ A% . 

4. Relation with directed-site lattice animals 

It was shown by Dhar in ref.0 that (d + 1) -dimensional problems of directed-site 
animals enumeration (DSAE) were related to <i-dimensional nearest neighbor exclusion 
models. The latter include the hard-dimer model (in d = 1 dimension) and the hard- 
hexagon model (in d = 2 dimensions), respectively related to DSAE on the simple square 
(d+1 = 2) and simple cubic (d+1 = 3) lattices. More precisely, the generating functions of 
DSAE were expressed as the density of occupation of the corresponding exclusion models, 
which in turn reduce precisely in the hard dimer and hard hexagon cases to the ratios 
yielding Z^°(t) and Z^°(t) (see eqns. ( |2.7| ) and (|3~^)). This strongly suggests the existence 
of a direct connection between semi-random Lorentzian lattices and directed-site lattice 
animals. In this section we will indeed establish an equivalence between the two problems. 
As an outcome, this will give an alternative derivation of the results of |§ in a more direct 
way. In the following, we will focus for simplicity on the cases of (1 + 1)-D Lorentzian 
triangulations and (1 + 2)-D Lorentzian tetrahedral complexes. 

4-1. From Lorentzian triangulations to square lattice directed animals 

Let us start with the set of configurations of vertical edges of an arbitrary Lorentzian 
triangulation, with a fixed left projection C made of hard dimers all occupying even time- 
slices. For each such configuration we first construct its skeleton, obtained as follows. We 
first decompose the edge configuration into blocks made of consecutive edges within the 
same time slice that are not separated by edges from neighboring slices (see Fig. 12 (a)). 
Squeezing each block into a single edge, we arrive at the skeleton of the configuration, itself 
a particular edge configuration with projection C, with no two consecutive edges within the 
same slice, as illustrated in Fig. 12 (b). Summing over all weighted configurations sharing 
the same skeleton simply amounts to assign an effective fugacity x = t/(l — t) for each edge 
of the skeleton. Indeed, x is nothing but the sum over all block sizes n = 1, 2, ... with a 
weight t n . We may now arrange the edges of the skeleton by representing its successive left 
projections along regularly spaced vertical lines numbered 1, 2, 3, ... (distant by ^/2, see 
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(a) (b) (c) 

Fig. 12: The equivalence between Lorentzian triangulations with a fixed left 
projection (with edges on even time slices) and directed site animals on the 
square lattice. We first build from the configuration its skeleton obtained by 
squeezing all blocks into single edges (a). We then place the remaining edges 
of the skeleton onto a square lattice tilted by 45° by pushing each edge to the 
leftmost available position on the lattice (b). This is equivalent to placing the 
successive left projections of the skeleton onto regularly spaced vertical lines 
here labeled 1,2,3,.... We have also indicated by arrows the two successors 
(in the next plane) of a vertex of the square lattice. A site on line % can 
be occupied by an edge only if one of its two predecessors (on line i — 1) is 
occupied. The occupied sites form a directed site animal on the tilted square 
lattice (c). 

Fig. 12 (b)). Noting that successive projections alternate between odd and even positions 
for all the edges, we see that they now lie on a regular square lattice (tilted by 45°). 
Moreover, the position of the edges satisfy the directed lattice animal constraint that a 
site on the vertical line i can be occupied only if one of its immediate neighbors on the 
lattice at vertical line i — 1 is occupied (see Fig. 12 (c)). In conclusion, the skeletons of 
configurations with projection C are in one-to-one correspondence with directed lattice 
animals on the square lattice with the same set C as "source" . We therefore end up with 
an identity between the generating function Aq{x) of square lattice animals with source 
C and activity x = t/(l — t) per occupied site and the partition function of Lorentzian 
triangulations with activity t per dual edge and fixed left projection C. 

A c (x) =t^ Z {c \t), witht=^— , (4.1) 

1 + x 

where we have used the thermodynamic partition function Z^ c \t) = liniT^oo (t), 
well-defined for any finite C. 
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The inversion relation (2"77) allows then to write Aq{x) as the density of occupation 



pc{z) of the set C in the hard dimer model with fugacity z = — t = —x/(l + x), namely 

yhd (_f\ 

Ac(x) = { -l)\c\ { - t) \c\^_L = (_i)|C| poH) (4.2) 



with t = x/(l + x) as above. This coincides with the result of Dhar 0. Our derivation 
using left projections is slightly more direct than that of @, which uses an intermediate 
connection with a dynamical crystal-growth model and relies on the similarity of master 
equations rather than on a direct correspondence between configurations. Note also that 
our inversion relation ( |2.5|) is not a simple rephrasing of the connection between directed 
animals and hard dimers, as it gives for instance access to the (finite size) partition function 
of the model of Lorentzian triangulations, which is not a density of occupation and has 
no direct animal interpretation. Moreover, in this section, we have artificially restricted 
our model to left projections with only edges in even positions, which was needed for the 
animal interpretation. However we have a Lorentzian surface interpretation for all other 
choices of the left projection C, and a relation to hard dimer densities too. 

4-2. From Lorentzian tetrahedral complexes to cubic lattice directed animals 

The above construction generalizes nicely to the case of Lorentzian tetrahedral com- 
plexes, which are in one-to-one correspondence with directed site lattice animals on the 
simple cubic lattice. We may repeat the above construction. We start again from a pla- 
quette configuration with fixed projection C made only of plaquettes lying on one of the 
three triangular sublattices @ numbered 1, 2, 3 of the triangular lattice formed by the cen- 
ters of the tubes. The skeleton of the plaquette configuration is obtained again by first 
decomposing it into blocks of successive plaquettes within the same tube, and shrinking 
each of these blocks into a single plaquette. Again, the partition function of the original 
configurations with fixed projection C and that of the skeletons with projection C are 
identified by assigning an effective fugacity x =t/(l — t) to the plaquettes of the skeletons. 
The plaquettes of the skeleton are then placed at the vertices of a simple cubic lattice by 
taking successive partial projections as follows. We start from the left projection C whose 
plaquettes lie on the sublattice, say 1, and place them on a first section plane numbered 1. 



The triangular lattice is naturally decomposed into three triangular sublattices, such that no 
two sites of the same sublattice are adjacent. 
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Fig. 13: The possible positions (a) for plaquettes in successive partial pro- 
jections of the skeletons. Plaquettes in planes 1, 2, 3, 4... lie in the sublattices 
1,2,3,1... The positions form a simple cubic lattice whose cell (b) is repre- 
sented together with the label of the corresponding [111] crystalline planes. 
The arrows in both pictures indicate the six immediate successors (in the 
planes 2 and 3) of a vertex of the plane 1. For a given skeleton, only part of 
the available positions are filled by plaquettes. A position of the cubic lattice 
can be occupied by a plaquette of a skeleton only if one of its six predecessors 
is occupied too. 

The left projection of the rest of the skeleton has plaquettes now lying on the sublattices 
2, 3 only. We pick the partial projection made of the plaquettes of this projection belonging 
to the sublattice 2, and place them on a second plane, numbered 2, parallel to number 1, 
and distant by l/y/2. Note that this partial projection might be empty. We then go on 
by considering the rest of the skeleton, whose left projection is now made of plaquettes 
on sublattices 3 or 1, out of which we retain those on the sublattice 3, etc ... Note that 
no two consecutive partial projection planes can be empty. The parallel partial projection 
planes 1,2,3... are nothing but the successive [111] crystalline planes of the simple cubic 
lattice (see Fig. 13 (a)-(b)). The skeleton may now be interpreted as a configuration where 
the vertices of the simple cubic lattice may be occupied (or not) by plaquettes, with the 
directed site animal constraint that a vertex in plane j may be occupied only if one of 
its preceding first or second nearest neighbors (respectively in planes j — 1 and j — 2) is 
occupied. This allows to identify the generating function for DSAE on the simple cubic 
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lattice with a fixed source C and an activity x per occupied site with the partition func- 
tion of (1 + 2)-D Lorentzian tetrahedral complexes with left projection C and a weight 
t = x/(l + x) per diamond-shaped building block (i.e. per tetrahedron). 

This in turn allows to recover the relation between DSAE on the simple cubic lattice 
and the hard hexagon model in two dimensions 0, in a slightly more direct manner. 



5. Transfer matrices and integrability 

In ref. Q it was shown that the existence of a transfer matrix formulation of pure 
Lorentzian triangulations was instrumental for deriving an effective one-dimensional 
(Calogero) Hamiltonian for the corresponding continuum scaling limit. It turns out that 
all the models introduced so far also have a more involved but similar transfer matrix 
formulation in terms of the original Lorentzian surfaces. In this section, we first show 
how to construct these transfer matrices, using in particular our inversion formula. We 
then show how these matrices can be regarded as particular points of integrable families 
of commuting matrices. 



5.1. Transfer matrices for Lorentzian surfaces 
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1 




-yVt 



-xvt 



(a) 



(b) 



Fig. 14: The transfer matrix of Lorentzian triangulations (a) connecting i 
lower half-edges to j upper ones across a given time-line. The corresponding 
(non-empty) hard-dimer configurations (b) and their respective weights. 



Let us first recall the form of the (infinite) transfer matrix for Lorentzian triangula- 
tions. This matrix T(t) transfers a row of half-edges across a time-line into another row 
of half-edges (see Fig. 14 (a)). A state is simply characterized by the number % (resp. j) 
of lower (resp. upper) half-edges, and the matrix element simply counts the number ( l ~^ J ) 
of ways of arranging these half-edges along the time-line, together with a weight yt per 
half-edge. Another way of obtaining this matrix T(t) consists in computing its generating 



29 



function @i(x,y\t) = 



i,j>0 



x 



y^T(t)i j by use of our inversion formula 



namely by 



expressing it as the inverse of the partition function of hard dimers on [0, 2], with weights 
—\ftx (resp. —y/ty) per dimer in the slice! (resp. 1) (see Fig. 14 (b)): 



Qi(x,y\t) 



1 



1 iy 



1 Wi 
iyft 



l-Vt(x + y) 



i,j>0 



t ]t — x l y J (5.1) 



Note that we have used a symmetrized version of the transfer matrix for hard dimers on 
the line, differing from that of ( 2.1 Q ) by a simple conjugation. The factors of i ensure that 
the dimers get the correct negative weight. This is appropriate to account for the fact that 
we count only half of the edge weight for the half-edges. With this definition, we have a 
similar formula for the T-th power of T(t) generated by @T(x,y\t), namely 

1 



e T (x,y\t) 



(i W. 



i \Vt 

iVt 



(5.2) 



corresponding to the inverse of the partition function of hard dimers on [0, T + 1] with 
a weight —t per dimer on the slices 1,2, ..,T — 1 and weights —\ftx (resp. —yty) in the 
slice (resp. T). This quantity is simply related to the loop-loop propagator ( [2.22| ) by 
e T (x,y\t) = G T +i(x/Vi, y/Vi\t). 




(a) 



(b) 



Fig. 15: The transfer matrix of Lorentzian surfaces made of time-like 
hexagons i.e. with dual edges of fixed length k = 2 (a) connecting lower 
and upper half-edges across a given time-line. The lower and upper half- 
edges each come in two types according to their position (1 or 2) along the 
full edge. The corresponding (non-empty) hard-trimer configurations (b) and 
their respective weights. 



We denote by m the slice [m, m + 1]. 
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Let us now turn to the model of Lorentzian surfaces made of 2(k + l)-gons defined in 
Section 2.3. To get a nice row-to-row transfer matrix Tk(t), we must distinguish between 
the possible states in which a given half-edge can be according to its relative positioning 
along the full edge of size k it belongs to (see Fig. 15 (a)). Let i r (resp. j r ) denote the total 
numbers of lower (resp. upper) half-edges occupying the position r = 1, 2, k along the 
edges of length k they belong to. The corresponding transfer matrix elements read 



2fe(i)i ll ...,i fc ;ii, 



,3 k 



k-l 

n*. 

r=l 



Jr-\-l i^r 



(5.3) 



expressing that a half-edge in position r transfers into one in position r+1 for r = 1, k— 1 
whereas one must choose the position of the j± new upper half-edges wrt all the other 
previously existing ones, with a weight per half-edge. Note that for k > 1 this transfer 
matrix is no longer symmetric. As before, we use the inversion formula ( |2.16| ) to write 
the free boundary partition function 6^ (x±, x k ; y±, yk\t) over a time lapse T as the 
inverse of the partition function of hard (k + l)-mers on [0, T + 1] possibly cut on the time 
extremities, with a weight —t per (k + l)-mer visiting the slices 1, 2, .., T — 1, and a weight 

2fc-2r+l 2r-l . . . 11/ m\ 

—x r t 2fe (resp. — y r t ^ ) for a [k + l)-mer cut in position r in the slice (resp. T), 
for T > k. When T < k, we must pay attention to (k + l)-mers that are cut at both 
extremities: these receive a weight —x r yT+ r t~ , r = 1,2, ...,k — T (see Fig. 15 (b)). This 
finally leads to 
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a 2k 1 Xk 



at ^k \ 






aH k 



(5.4) 



/ 



where we have set a = e«. The states on which Tk(t) acts can be understood as being 
respectively the vacuum, the first (lowest) elementary segment of the (k + l)-mer, the sec- 
ond, third,..., k-th one. As before, the phases in ( |5.4| ) ensure that all multimers (including 
the cut ones) have a minus sign in front of the fugacity. Note also the explicit "up-down" 
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symmetry of Tk(f) implemented by the symmetric matrix IZij = dj^+i-i for i, j > 1 and 

TZqj = 8j t o, namely that the matrix 1ZTk(t) is symmetric and that IZvk(x) = Wk(x). This 

(k) 

translates into the symmetry of 0^ under the interchange of x r <-> y^+i-r, r = 1, 2, k. 



A preliminary remark is in order: if we take the formula ( |5.4| ) for T = 0, we get 



8? ) (xi,...,x fc ;y 1 ,...,y fc |t) = — — jj = £ fl^)" ( 5 " 5 ) 

which generates the diagonal matrix P with diagonal entries Pi lt ... y i k = ^ ! rather than 
just the identity. More generally, ^ is the generating function for (Tfc(t)) T P, where the 
boundary operator P corresponds to a preliminary ordering of the half-edges of the initial 
state over which Tk(t) acts. Note that the action of Tk(t) automatically orders the final 
state. In particular, for T = 1 the result ( |5.4j ) reduces to 



9^(aPi,...,aPfe;i/i,...,i/jfe|t) 



l-t^k ( Xfc + yi) -ti(x x y 2 + x 2 y 3 + ■■■ + Xk-iVk) 

h,... ,ifc,ii>0 

which is nothing but the generating function of Tk(t)P, with Tk(t) as in ( p.3| ). The up-down 
symmetry is recovered by noting that RTk(t)R = P~ 1 Tf-(t) t P, where R is the matrix with 
entries Ri 1 ,...,i k -,j 1 ,...,j k = H r= i S jr>ik+1 _ r corresponding to the interchange x r <-> x k +i- r , 
while the transposition expresses the interchange of x r and y r . 

More generally, the above construction can be generalized to the hard multimer models 
of Section 2.2, leading to an explicit transfer matrix formulation of the corresponding 
generalized Lorentzian surfaces. Details are gathered in Appendix A below. Finally, let 
us mention that the above extends straightforwardly to higher dimensions. In the case 
of (1 + 2)-D Lorentzian tetrahedral complexes for instance, we may define a "plane-to- 
plane" transfer matrix acting in the second time direction (T2, see Fig. 10) and transferring 
from a (1 + 1)-D Lorentzian triangulation to the one above (see Fig. 11). The generating 
function for this transfer matrix (multiplied by a diagonal ordering operator P as before) 
is expressed as the inverse of a decorated partition function for hard hexagons via our 
inversion formula. The latter is itself the generating function for the matrix elements of 
the row-to-row transfer matrix for hard hexagons 0]. 
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5.2. Integrability and the inversion relation 

In ref. it was observed that the model of pure Lorentzian triangulations was a point 
(a = 1) along an integrable family of models including both triangles and squares (with 
two time-like and two space-like edges), with weights ay/i per triangle and u = til — a 2 ) 
per square. The corresponding surfaces are generated in the dual picture by means of 
a row-to-row transfer matrix T(t, a) whose element T(t, a)ij corresponds to the transfer 
from i lower vertical edges to j upper ones across a time-line, i,j > 0. It was shown that 
the transfer matrices T(t, a) and T(t', a') commute provided 

l-t(l-a 2 ) = l-t'(l-a /2 ) (5 7) 

aVt a'Vt 7 

In our language this model with triangles and squares corresponds to having 2(k + l)-gons 
of arbitrary size k = 1, 2, 3, ... made of the piling up of {k — 1) squares across k consecutive 
time-slices, and terminated by two triangles on the top and bottom of the pile. These 
objects come with a weight tk = a 2 tu h ~ 1 , k = 1,2,3,.... We immediately deduce from 
our general inversion formula that the partition function fi(t, a) for these surfaces with 
left staircase boundary conditions and with semi-infinite time range [0, oo) satisfies the 
equation ( |2.20|) 

oo 

fi = l + 5]aV"V H1 t^ 2 = (l + u)^t-l (5.8) 

k=i 

Upon the rescaling fl = (1 + u)fi and t = t/(l + u) 2 we end up with the same equa- 
tion as that of the pure Lorentzian triangulation case a = 1 (u = 0). As a non- 
trivial outcome of our formulation, the integrability condition for the infinite size transfer 
matrices T(t, a) can be deduced from a similar condition on the 2x2 transfer matri- 
ces of the corresponding hard multimer models. Let us write the generating function 
02(x, y\a, t, a', t') = J2i j>o xl yHT(t', a')T{t, a))ij of the product T(t' , a')T(t, a) as the in- 
verse of a multimer partition function on [0,2] with weights —a\/tx (resp. —aa'Vtt' and 
—a'y/Fy) for dimers in slice (resp. 1 and 2), —a'Vt/ux (resp. —ay/tu'y) for trimers cov- 
ering the slices 0,1 (resp. 1,2) and finally —xyuu' for quadrimers covering all the slices. 
This reads 

6 2 (x, y\a, t, a', t') = — J. : — (5.9) 

, A . , I 1 la yt' \ I 1 la^/t \ I 1 

( 1 iy 



10, Vt 7 u' ) \ iaVi 
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Similarly the generating function 2 (x, y|a', t', a, t) is expressed by the same formula with 
the two 2x2 matrices exchanged. The commutation of the original transfer matrices 
translates therefore into that of these two 2x2 ones. A simple calculation shows that they 
commute provided (1 — u)/{a\fi) = (1 — u')/(a'\rF), which is precisely the integrability 
condition (|5.7|). 

This suggests to look for integrable families of finite size matrices corresponding to 
decorated multimer models and to interpret them as integrable models of more involved 
Lorentzian surfaces. For simplicity, let us consider the case of symmetric (k + 1) x (k + 1) 
matrices. For illustration, we may consider a model of colored hard dimers with a color 
index i = 1,2, ...,k, and a weight z% for each dimer of color i, with transfer matrix 



z{ ••• 



zi ■•• 



(5.10) 



■•• / 

This is nothing but the original hard dimer model with fugacity z = z\ + z<i + ... + z^. 
However, we now want to view it as a particular point of an integrable family of transfer 
matrices T, depending on k(k + 3)/2 new parameters Tj = Tj^, i,j = 0, 1, 2, k where 
we fix Tq q = 1. The parameters %j = y/wj, j = 1, 2, k are interpreted as new fugacities 
per dimer of color j. The numbers %j for i,j>l can be interpreted as follows. We first 
allow the dimers to pile up so as to form multicolored multimers of length m, and we assign 
a weight ^/Wi x Wi m EI^Li 1 % v ,i v +i f° r each (m + l)-mer formed by a sequence of dimers of 
colors ...,% m - To get an integrable family, we must impose that all the T share the 

same orthonormal basis of eigenvectors as Tq (this involves fixing a (k — 1) -dimensional 
rotation acting on the kernel of Tq) . Denoting by the z-th entry of the m-th eigenvector 
of To, and by A m the corresponding m-th eigenvalue of T, for i, m = 0, 1, 2, k, we have 

T hJ = £ A m 4 m) ^ m) (5-11) 

m=0 

where the A m 's are free parameters satisfying the constraint 

r ,o = l = ^A m (4 m) ) 2 (5.12) 

m 

This shows that the matrices T form generically a /c-parameter family passing by the point 
Tq (with Wi = Zi). This also gives a parametrization of the dimer factors Tq^ = ypw[ = 
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J2 m ^m^^ "ipt"^ m terms of the k + 1 eigenvalues A m . Let us display the case k = 1 
for illustration. We first diagonalize the matrix Tq = (^fj. *"q^ with the resulting 



eigenvectors: 



^ ' = _L= (1, i g ) , 4 = fa -1) (5-13) 



where we have set 1/y/t = q + 1/q, and q < 1. The one-parameter integrable family takes 
the form 

Tij = (1 + i ?v ^)^ ( °V 7 (0) + (1 + iy/w/q)^ 1 ^ (5-14) 



where we have parametrized the eigenvalues Aq and Aj, satisfying Ao(V ; o°^) 2 +Ai( , i/>Q 1 ' ) ) 2 = 1 
in terms of the matrix element Tqi = y/w. Eqn.( j5?14|) coincides with T = ^ io-V^ \ 
provided = iay 7 ^ and w = £(1 — a 2 ). 



A slight generalization of the previous case allows to include the general Lorentzian 
surfaces models of Sections 2.2 and 2.3 in larger integrable families. As opposed to the 
previous case, the corresponding transfer matrix Tq is no longer symmetric (see for instance 
eqn. (|5.4|) ). Still, the "up-down" symmetry of the problem allows to make the transfer 
matrix symmetric, by use of a certain involution i — > r{i) of its indices, implementing 
the up-down reflection of the multimers, while preserving the vacuum (r(0) = 0). More 
precisely, let 1Z be the matrix with entries IZij = Sj >r u\, then the up-down symmetry of Tq 
translates into the fact that 1Z Tq is symmetric. We now may look for an up-down symmetric 
integrable deformation T of Tq, namely such that JZT is symmetric, and moreover satisfying 
the condition that 7q,o = 1- Denoting again by a basis of diagonalization of Tq, we 
first note that by virtue of up-down symmetry it can be normalized so as to be U 1Z- 
orthonormal" , namely that (TZi/j^Y i/^^ = 8 m , p . We still write the family T in a form 
similar to ( |5.11| ), taking into account the up-down symmetry, namely 

^•=]TA m ^V^ (5-15) 

m 

where the A's satisfy the same constraint as before: Tqq = 1 = J2 m ^mhPo"^) 2 ■ The case 
of hard (k + l)-mers of Section 2.3 is obtained by taking To = Tk(t) of eqn. (|5.4|), with 
the same up-down symmetry matrix TZij = Sj^+i-i for i,j = 1,2, k and TZqj = 8j t Q. 
The case of hard multimers of Appendix A corresponds to taking Tq = T^ hm \ and the 
up-down symmetry matrix 1Z = TZ( hrn ^ . 
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In conclusion, we have shown that all the ID hard-multimer transfer matrices used in 
this paper are part of larger integrable families corresponding to the introduction of extra 
interactions between multimers. This in turn implies that all the transfer matrices for 
Lorentzian surfaces defined in Section 2 are also part of larger integrable families allowing 
for the original polygonal tiles to form larger composite objects, with new weights. This 
extends nicely to the case of (1 + 2)-D Lorentzian tetrahedral complexes as well. Indeed, 
as noted by Baxter [[7], the hard hexagon model is part of a larger integrable family (hard 
square model with diagonal interactions). Once translated in the language of plaquettes 
this amounts to allow for hexagons to coexist on neighboring sites along one of the three 
directions of the triangular lattice of tubes, with a new contact Boltzmann weight. The 
corresponding transfer matrices form via our inversion relation the desired integrable family 
of models, describing complexes in which the diamond-shaped building blocks can be glued 
together along that same direction so as to form larger straight objects of arbitrary length, 
with new Boltzmann weights. 

6. Discussion and conclusion 

In this paper, we have first investigated general models of (1 + 1)-D Lorentzian surfaces. 
Viewing the pure Lorentzian triangulations as semi-random surfaces made only of time-like 
lozenges, we have considered semi-random surfaces made of larger time-like tiles, namely 
obtained by piling up lozenges, thus forming 2(i + l)-gons, i = 1,2, 3.... All these models 
have been solved by use of an inversion relation expressing their partition functions in 
terms of those of associated hard multimers in ID. We have then showed how to obtain 
multicritical models of these generalized Lorentzian surfaces by an appropriate fine-tuning 
of the activities of the various tiles. These models involve only time-like interactions, 
a key point for the inversion relation to ID hard objects to hold. Indeed, this relation 
relies entirely on the fact that the only horizontal (space-like) coupling is through the no- 
crossing prescription of edges. It would be interesting however to try to build more involved 
inversion relations, possibly using horizontal (space direction) transfer matrices between 
successive vertical projections, and use them to solve models with horizontal (space-like) 
as well as vertical (time-like) interactions, such as the Ising model on Lorentzian surfaces. 

We then generalized our construction so as to build a new higher-dimensional model of 
(1 + 2)-D Lorentzian tetrahedral complexes, with two regular (time-like) and one random 
(space-like) directions. This model was solved by using a generalized higher-dimensional 
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version of the previous inversion formula, relating it to the 2D hard hexagon model solved 
by Baxter ]7[]. This model is the most natural extension of that of Lorentzian triangula- 
tions to higher dimension with tetrahedra playing the role of triangles. We also showed 
how to interpret the model as the time evolution of Lorentzian triangulations in successive 
layers. In this respect, we obtain a toy-model for (2 + 1)-D Lorentzian gravity, where the 
space-like triangulations are restricted to be themselves Lorentzian. It is interesting to 
note that other plaquette models can be similarly constructed by considering tubes with 
sections forming different lattices. For instance, we may consider a model of square plaque- 
ttes with section the square lattice or of triangular plaquettes with section the triangular 
lattice. In an analogous way, the latter are respectively related to the hard square and 
hard triangle models in two dimensions. Unfortunately, these have not yet been solved 



exactly, although very precise results are known about them [[16| . More importantly, when 
we try to interpret these as models of simplicial Lorentzian manifolds, we find that the 
corresponding elementary simplices (replacing the tetrahedra) become degenerate objects 
with pairs of vertices linked by more than one edge, that are usually discarded in simplicial 
gravity. This is why we did not go into details of these models. Another possible gen- 
eralization is to keep hexagonal tubes and tetrahedra as elementary objects, but to form 
larger building blocks than the diamond-shaped dodecahedra considered so far. These 
larger objects are obtained by gluing several dodecahedra in the two time directions in 
the same way as we did for lozenges. Models based on these building blocks with say 
activities ti for the i-th type of block are still related to models of hard objects on the 
triangular lattice, namely with plaquettes obtained by gluing hexagonal plaquettes to- 
gether, and attaching to them activities Z{ = —ti. These 2D models seem not to have 
been considered yet (in particular, they are not identical to the so-called ABF models 
|17|| generalizing the hard hexagon model). Still it is very plausible that by again fine- 
tuning the activities ti one should be able to reach multicritical points. It is natural to 
believe that these correspond to non-unitary conformal field theories with central charges 
c(2, 2m — 1) = 1 — 3(2m — 3) 2 /(2m — 1), where the case of hard hexagons corresponds 
to the Lee- Yang edge singularity with m = 3, c = —22/5. The corresponding thermal 
dimension is known to be = (5 — 2m) /(2m — 1) henceforth the singularity of the 
free energy has the critical exponent 2 — a = 1/(1 — ^1,3) = (2m — l)/(2(2m — 3)). The 
corresponding fractal dimension of the associated multicritical Lorentzian manifolds would 
read d F = 4(2m - 3)/(2m - 1). 
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Let us now point out that there is an obvious natural higher- dimensional generalization 
of tetrahedral complexes, now built with d+ 1-simplices living in d+ 1-dimensional "tubes" 
whose (i-dimensional section is the regular (i-dimensional generalization of the FCC lattice 
||18j| , in correspondence with a model of nearest neighbor exclusion on the vertices of that 
lattice. 

Another direction of generalization consists in staying in low dimension e.g. in (1 + 1)- 
D, but to introduce disorder in the form of randomly distributed activities, but constant 
within each tube. Indeed, the inversion relation of this paper holds as well for each realiza- 
tion of the disorder, with fixed tube-dependent activities. In the case of quenched disorder, 
where we must compute the average Log Z over all possible realizations of the disorder, 
we can use the inversion formula term by term in the average, to write the result as minus 
the quenched disorder average of the corresponding hard-object model, namely — LogZ^. 
Clearly, the inversion relation does not allow to relate the two annealed disorder partition 
functions Z and Zh of the two problems. 

Another type of disorder consists in considering plaquette models with fixed activities 
but defined on tubes whose section is itself arranged into random graphs, in connection 
to the corresponding nearest neighbor exclusion problem on the same graphs. Again, the 
two quenched disorder problems are identical up to a sign. 

We may also consider tubes with section arranged into an arbitrary but fixed graph. 
An interesting choice for such a graph is one where the nearest neighbor exclusion model 
has already been solved, in which case we can immediately convert the solution into that 
of the corresponding Lorentzian plaquette model. One such example is the infinite rooted 
g-valent graph known as the Bethe lattice, q = 2, 3, .... The corresponding plaquettes are 
g-gons, each connected to one q-gon in the previous shell and to q — 1 ones in the next. A 
simple calculation (see e.g. ||16|| ) shows that the thermodynamic free energy of the nearest 
neighbor exclusion model f n ne{z) per site, with an activity z per occupied site, reads 

f nne (z) = (1 - |) Log(2 - ijl(-z)) - | Log/i(-z) (6-1) 

where //(£) is the Fuss-Catalan generating function ( |2.24| ) with k = q, obeying the equation 
z = (1 — n(t)) / (i(t) q , t = —z. For q > 2, the first singularity of f nn e(z) for z < occurs 
at the value z = z c such that z c ) = q/(q — 1), resulting in a singularity with exponent 
2 — a = 1/2, translating into a fractal dimension dp = 2 for the corresponding Bethe 
Lorentzian surfaces. 
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Finally, the models considered here were shown to contain lattice animals as a subclass. 
Conversely, we may apply this equivalence backwards to investigate more refined properties 
of our semi-random lattices. In particular, directed lattice animals are known to have two 
distinct characteristic lengths, a longitudinal and a transversal one, each with its own 
scaling exponent v\\ and u±. While u± refers to our time correlation length exponent 
(with (If = l/v±), it would be interesting to interpret the second exponent in terms of 
our semi-random lattices. On the other hand, the animals interpretation holds as well for 
models with 2(i + l)-gons % = 1, 2, 3.... However the animals corresponding for instance to 
(1 + 1)-D models with edges of fixed length k > 2 are quite unconventional, not to speak 
about those corresponding to the multicritical models. 

After completion of this work, we became aware of ref. []T!| where an analogous in- 
version formula as that for (1 + 1)-D pure Lorentzian triangulations was derived in the 
context of commutative monoids, and of more recent developments on lattice animals and 



exclusion models []2(| [pT 
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Appendix A. Transfer matrix for general multimers 

Let us consider the general model of Lorentzian surfaces made of 2{i + l)-gons i = 
1, 2, k of Section 2.2. It is easy to see that the corresponding transfer matrix T^ hrn ^ is 
expressed in terms of a hard multimer one (denoted by T^ hm ^) through a formula similar 
to ( |5.4|) . The matrix T^ hm ^ is obtained by superposing those for fixed i (with respective 
transfer matrices Ti(t) as in (|5.4j) ) in the following manner. It must act on the following 
states: the vacuum (label 0), the unique dimer (label 1), the first segment of a trimer 
(label 2), the second segment of a trimer (label 3),..., the j-th segment of a (i + l)-mer 
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(label = j + i(i — l)/2 for j = 1, 2, i),... the k-th segment of a (k + l)-mer (label 

/c(/c + l)/2). The matrix 7~(^ m ) has therefore total size 1 + k(k + 1)/2 and its elements read 



rr-ihm) _ 1 

-*0,0 — L 

^0 = ^(^,0 *.J>1 (A- 1 ) 
^I(i,j),I(i,m) = Ti{t)j tm i,j,m > 1 
/^-(/i?n) — — q for 2 ?5 

where Tj(£) is as in ( |5.4j) . Note that this matrix is also "up-down" symmetric in the sense 
that Ji( hm )']'( hrn ) j s symmetric, where TZ^ hm ^> is the symmetric matrix with entries 

•nW _ i 
71 -o,o — L 

Kffi .s = TZ^X n = (A.2) 



<rp(hm) _ r j- 

'^I(i,j),I(p,m) ~ u p^ u m,z+l-j 

for 1 < j < i < k. The matrix T^ hrn ^ leads through a formula analogous to (15.4!) to 
the ge^attag taction B^iM feWk. U) of (T<->fp<-> whtfthe 
diagonal matrix p( hm ) implements the ordering of the initial state over which T^ m ) acts. 
More precisely, we have 



T ({^/(i,j)};{2//(p,m)}|^l, -,**) 



(A 3) 

with wr(;j) = (wi)j(yi(i,j)), w = l 
vi(i,j) = {y i ) j {x I{i j ) ), v = i 

which for T = generates the elements of the diagonal matrix p( hm ) with diagonal entries: 

T-)(hm) 

^(,,-)} n^( P ,m)! 1 ' 1 

while the formula ( |A.3[ ) leads for T = 1 to the transfer matrix elements 

k / p-1 

{vcp.m)};{j/(p,™)} 111 ll°ji(p, r +i),»j(p, r ) 

p=l \ r=l 



(A.5) 

(Ep=l(j/(p,l) +Er=l'%))) ! 
(E l<r<p<A; Ilp=l i/(p,l)' 
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For illustration, the matrix T^ hrr ^> reads for k = 2: 



/ 1 















i 








V 








rj-(hm) 



where a = e 17r//4 . This leads to the matrices 

p(ta) (h +i2 + hV- 



«*2 \ 





/ 



H ,*2,*3 



-,{hm) 



"»l,«2,«3;jl,j2,j3 



2 '2+ i 3+J2 

U 2 U 4 5 



33^2 



[h +12 + 13 + jl +J2) 

(h +i2 + hV-ji } -j2 } - 



(A.6) 



(A.7) 



where i\ (resp. £2,^3) denote the total numbers of halves of single edges (resp. halves of 
first, second segments of edges of length 2) in the lower part and similarly for the j's in the 
upper part. Note that the combinatorial factor in the second line of ( |A.7|) expresses the 
choice of position of the two new types of edges (of respective length 1 and 2, in numbers 
ji and j'2) wrt the already existing ones. 
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